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Abstract. A wide class of boundary problems in quantum mechanics is 
discussed by using path integrals. This includes motion in half-spaces, radial 
boxes, rings, and moving boundaries. As a preparation the formalism for 
the incorporation of 5-function perturbations is outlined, which includes the 
discussion of multiple 5-function perturbations, 5-function perturbations along 
perpendicular lines and planes, and moving 5-function perturbations. The 
limiting process, where the strength of the 5-function perturbations gets infinite 
repulsive, has the effect of producing impenetrable walls at the locations of the 
5-function perturbations, i.e. a consistent description for boundary problems 
with Dirichlet boundary-condition emerges. Several examples illustrate the 
formalism. 



1. Introduction 

Boundary problems, respectively boundary- value problems, appear at almost every 
stage of quantum physics. In a conventional treatment boundary-conditions are 
usually used to define the solution of a particular solution, say of the Schrodinger 
equation, in an unambiguous way. Such boundary-conditions, e.g. vanishing wave- 
functions at infinity, are the basic requirement for the very set-up of the relevant 
Hilbert space. 

Boundary-conditions at infinity, however, are a very specific idealization and 
for many more physical situation not appropriate. Typical experimental situations 
require boundary-conditions at a finite distance from the origin, for instance in 
electrodynamics Dirichlet or Neumann boundary-conditions are the most discussed. 
Analogous considerations in quantum mechanics lead to Dirichlet and Neumann 
boundary-conditions, respectively, for the wave functions. In the following we will 
consider such quantum mechanical problems with Dirichlet boundary-conditions. 

To deal with these problems, there have been only few attempts in the literature. 
For instance, numerical studies can be found in Refs. [2, 14, 54, 62] for the problem 
of an oscillator in a box. 

Studies of magnetic systems in euclidean and radial boxes in the context a solid 
state physics are due to Falkovsky and Klama [22, 23], and recently by Klama and 
Rossler [46] . In these studies the authors succeeded in the explicit construction of the 
relevant Green function by exploiting the given boundary-conditions for a particular 
system. However, their approach lacks a systematic approach to the whole problem 
of boundary-conditions, in particular for multidimensional boundary problems with 
several boundaries. Also, path integral considerations were not made. 

Path integral studies for boundary problems have been done by several authors. 
First, by the explicit construction of the propagator by using the "mirror" principle 
(see below), c.f. in the classical textbooks on path integrals of Feynman and Hibbs 
[25] and Schulman [57], and more recently by Inomata and Singh [44], and Janke and 
Kleinert [45], and second by an attempt to build in general boundary conditions into 
the path integral, see e.g. Clark et al. [13], and Carreau et al. [9, 10]. The best known 
examples are that of bounded, but otherwise free motion, i.e. motion in half-spaces, 
euclidean and radial potential wells, which are all standard examples. 

The explicit construction of the propagator for boundary problems with, say, 
Dirichlet boundary-conditions is generally quite involved. For simple systems, how- 
ever, the mirror "principle" can be applied. Let us consider the boundary at x = as 
totally reflecting, and denote the reflecting action by the operator 7. The propagator 
corresponding to a quantum Hamiltonian H (x G M) 




(1.1) 
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of the time-evolution equation 

V(x",t")= [ K (x",x';t",t , )^(x , ,t , )dx' (1.2) 

then can be used to construct the propagator in the half-space IR + in the following 

way: 

K* i+ (x",x';T) = K (x",x';T)-K (x", 1 x';T) . (1.3) 

Here it must be provided that Kq(T) has the invariance property Kq^x", 7a/ ; T) = 
Kq(x" , x';T), respectively the corresponding Hamiltonian is invariant under the action 
of 7: 7-ffo = Hq. Examples, in the case of the half-space IR + , are the free particle, 
the harmonic oscillator, or the potential V(x) = Vq/ cosh 2 x [32, 34]. 

For a (free) particle in a box, there is an infinite number of reflections and 
translations which must be summed over, which gives the propagator for a particle 
in a box in terms of Jacobi 6-functions, a well-known result [25, 57]. That is, the 
propagator K{T) is then constructed as 

K(x",x';T) = Y J K Q {x", 1 x';T) , (1.4) 

and T denotes the set of all actions 7. This general "mirror" principle which also 
includes group actions plays an important role in the theory of semiclassical periodic 
orbit theory and quantum chaos (see e.g. [36, 58, 60] and references therein). 

However, the conceptual simplicity of Eq. (1.4) may be striking, but in most cases, 
K(T) cannot be explicitly evaluated, or even Eq. (1.4) cannot be applied because H 
and Kq(T) are not invariant with respect to 7. 

In this paper I want to consider quantum mechanical problems with boundary- 
conditions for a wide class of systems. I will consider only Dirichlet boundary- 
conditions and do not discuss generalizations of more general local boundary-conditions 
along the lines of Refs. [9, 10, 13]. The idea starts by first considering a 5- function 
perturbation in the path integral [34] . The path integral can be expanded into a per- 
turbation series, the 5 functions allow to integrate over all intermediate coordinate 
positions, and the convolution theorem of the Fourier transformation decouples the 
time-convolutions, giving a geometric power series in the perturbation which can be 
exactly summed, yielding finally a closed expression for the (energy-dependent) Green 
function of the perturbed problem (c.f. the monograph by Albeverio et al. [1], as the 
main reference for a sound mathematical definition of such 5-function perturbations 
in quantum mechanics). 

Of course, by repeating this procedure, an arbitrary number of 5-function pertur- 
bations can be taken into account. 
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Dirichlet boundary-conditions at the location of the 5-function perturbation are 
now generated by making the strength of the perturbation infinitely repulsive [13, 52], 
thus yielding a closed expression for the Green function of the boundary problem. 
Having introduced the first boundary, we obtain by repeating this procedure with a 
second 5-function perturbation the Green function in a box with Dirichlet boundary- 
conditions on both boundaries. 

As I will show, it is possible to generalize this formalism to higher dimensions. 
Here an appropriate combination of 5-function perturbations along perpendicular 
lines, planes and hyperplanes can be introduced, which produce by the same limiting 
procedure infinitely repulsive walls at the location of the 5-function perturbations. It 
is obvious that arbitrarily shaped boundaries cannot be taken into account by our 
formalism. 

The further contents of this paper is as follows. 

In the next section I shortly sketch the theory of incorporating 5-function pertur- 
bations in the path integral. The whole theory was already presented in some length in 
Ref. [34] and is not repeated here in all details. However, I will develop the formalism 
further and will discuss 

1) Multiple 5-function perturbations. Closed formulas for the Green function for such 
arrangements were were not presented in [34] (compare also Ref. [24], however 
with the restriction to the wave-functions only), 

2) 5-function perturbations along perpendicular lines and planes, and 

3) Moving 5-function perturbations. 

These results are used in the third section to construct by a limiting procedure 
Green functions for boundary problems with Dirichlet boundary-conditions. I consider 
half-lines, one-dimensional boxes, half-spaces and two-dimensional boxes (with obvious 
generalizations to higher dimensions), and moving boundaries. 

In the fourth section I present several examples to illustrate the formalism. The 
examples will include 

1) Boundary problems in half-spaces, 

2) Boundary problems in boxes, 

3) Boundary problems in radial boxes, 

4) Boundary problems in radial rings, and 

5) Moving Boundaries. 

Numerous examples for 5-function perturbations in the path integral were already 
given in Ref. [34] and are not repeated here. 

The fifth section is devoted to a summary and discussion of the results; in Appendix 
A, the derivation of the Green function for multiple 5-function perturbations is given, 
and in Appendix B the Green function for the linear potential is constructed. 
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2. Path Integrals with 5-Function Perturbations 



2.1. Summation of the Perturbation Expansion 

The general method for the time-ordered perturbation expansion is quite simple. We 
assume that we have a potential W{x) = V{x) + V(x) in the path integral, where it is 
assumed that W is so complicated that a direct path integration is not possible. 
However, the path integral corresponding to V(x) is assumed to be known. We 
expand the path integral containing V(x) in a perturbation expansion about V(x) 
in the following way. The initial kernel corresponding to V propagates in At-time 
unperturbed, then it is interacting with V, propagates again in another At-time 
unperturbed, a.s.o, up to the final state. This gives the series expansion [25] (x G IR) 

x(t")=x' 



+ " r- 



K(x",x';T) = J Vx(t)exp^J yi 2 - V(x) - V(x) dtj 



OO s 

= K( v \x",x';T) + J2[ 




n=l 

xK { y\x u x';t l -t')V{x 1 )K { y\x 2 ,x l -t 2 -ti) x .. 
•••x V(x n -i)K {y \x n ,x n -i;t n -t„_i) x ... 
■■■xV(x n )K( v \x",x n ;t" -t n ) 

xK { y\x u x';t 1 -t')V{x 1 )K { y\x 2 ,x 1 -t 2 -ti) x ... 

• • • x V {x n -\)K {y \x n ,x n -\\t n - t n -i) x ... 

■■■xV(x n )K^ v \x",x n -t" -t n ) . (2.1) 

In the second step I have ordered the time as t' = to < t\ < t 2 < • • • < t n+ i = t" and 
paid attention to the fact that K(tj — tj-i) is different from zero only if tj > tj-i, 
([34], see also e.g. Bauch [6], Lawande and Bhagwat [50]), and we set T = t" — t' . 

We consider now an arbitrary potential V(x) in one dimension with an additional 
5-function perturbation so that [34] 

W{x) = V{x)- 1 8{x-a) . (2.2) 
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The path integral for this potential problem reads 



<t")= 



K^ 5 \x",x';T) 



Vx(t) exp 



m 



x — V(x) + r )8{x — a) 



dt 



(2.3) 



x{t')=x' 



We have assumed that the path integral (Feynman kernel, respectively) for the 
potential V is known, i.e. we set 



x(t")=x" 



Vx(t) exp ^ — 



rn 



K^ v \x",x';T) = 

x(t>)=x> 

including, of course, the (energy-dependent) Green function 



■-x 2 - V(x) 



dt 



(2.4) 



G^ v \x\x'-E)= l - / dTe iET / h K^(x",x';T) 
* Jo 

K^ v \x",x';T) = ^ e - iET / n G^ v \x'\x';E) 

J-oo 271-1 



(2.5) 



Introducing the Green function G^ S '(E) of the perturbed system similarly to (2.5), 
it is easy to sum up the emerging geometric power series, and we obtain due to the 
convolution theorem of the Fourier transformation 



G^(x",x';E) = G^(x",x';E) 



G ( y\x",a-E)G i y\a,x';E) 
G<Y)(a,a;E) -l/>y 



(2.6) 



where it is assumed that G^ v \a,a;E) actually exists. The energy levels E n of the 
perturbed problem W(x) are therefore determined in a unique way by the equation 



G^ v \a,a;E n ) = - . 

7 

The wave-functions are given by the residua of G(E) at E = E n , i.e. 



(2.7) 



* n (x) = lim 

±L — >£/„ 



E n — E 



l/^-G( v \a,a;E) 



1/2 



G {v \x,a;E) . 



(2.8) 



Here and in the following the index n denotes the level index of a corresponding energy 
level E n , with n — 0,1, ... , Nm, and Nm may be finite or infinite. 

Let us note that an implicit equation for the time-dependent propagator is due to 
Gaveau and Schulman [27]. They obtained 

t" 

K^\x",x';T) = K { y\x",x';T) + i^ [ K^ v \x" ' ,a;t)K^ 5 \a,x';T - t)dt . (2.9) 

« Jv 



Similarly as for one-dimensional problems, we can consider in D dimensions a 
radial potential according to 



W(r) = V(r) — 7<5(r — a) , 



(2.10) 



i.e. a spherical shaped 5-function (also called "shell"-, "surface"- ^-interaction [1, 
3]) perturbation of the potential V(r). We must require a ^ because a point 
perturbation leads to the evaluation of Green functions, where both arguments are 
equal (and zero), an expression which in general does not exist. Of course, we are 
using the usual D-dimensional polar coordinates. For a radial problem we can separate 
variables in the path integral: 



x(t )=x 



K^ w \x",x';T) 



x(t')=x' 



Vx(t) exp 



k Li 



rn 



x 2 - W{r) 



dt 



= "£Kr\r",r';T)S?(n")Sr(V) , 



(2.11) 



1=0 



where S^(Q) are the real hyper-spherical harmonics of degree / with unit- vector O [20] 
and the radial path integral K^ W ^ (T) is given by 



K\ W \r"y-T) 



r(t")=r" 



(r'r") 1 ^ j Vr(t)exp^£ 

r(t')=r' 
r(t")=r" 

(r'r") 1 -^ J n u n-2[r 2 ]Vr(t)exy 

r(t')=r' 



rn 



r 2 -h 



2 l(l + D-2) 
2mr 2 



- W(r) 



fl Li 



rn 



r 2 - W(r) 



dt 



dt 



(2.12) 



The functional weight \i v \r 2 \ is defined by by [39, 61]: 



N 



ri u [r 2 ] = lim ( '[n„[r j - 1 r j 



TV 



rT / 2Tvmr J - 1 r j \ 1/2 / mr J - 1 r j \ ( mrj^rj 



(2.13a) 



(2.13b) 



in order to guarantee a well-behaved boundary behaviour at initial at final points. I v 
describes a modified Bessel function. Proceeding similar to the previous section and 
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expanding the perturbed problem into a perturbation series yields: 



r{t")=r" 

^jf dTe iET/h J ^ + i^[r 2 ]Pr(t)exp|^ ^x 2 - V (r) + 7 5(r - a) dtj 

r(t')=r' 
oo 

= ^G{ 5) (r // ,r / ; J B) ( 5f(fi // )-5f*(n / ) , (2.14) 
where the radial Green function is given by 

Gf\r"^E) = Gl V \r",r';E) - G\ v \^E)G^( a y ; E) 

G ( l V \a,a;E)-l/a D - 1 <y 

Therefore the energy-levels E n are determined by the equation 



-J- T - = G ( l V) (a,a;E n ) . (2.16) 

The corresponding wave-functions are given in a similar way as in the one-dimensional 
case. 



2.2. Multiple 5 -Function Perturbations 

Starting from Eq. (2.6) we first consider the problem of two 5-function perturbations 
(such problems play e.g. a role in the study of the Casimir effect between to perfectly 
conducting planes [7]) 

K^\x",x';T) 

x{t")=x" 



x(t')=x' 



—x - V(x) + 7i 5 (a; - ai) + 72<5(x - a 2 ) 



(2.17) 



We obtain 



G<»V, x'; g) = G^V fl) - G< *'^" °* f E ) , (2 . 18) 



G( 5 i)(a2,a 2 ;E) 

72 



where G^ Sl '(E) in turn is given by 



C<* V, ^ E) = G<y\x", x'; E) - G(V){X "> ai; glg^W^ g) . (2 .i 9) 

G(^)(ai,a i; E)-- 
7i 



Insertion of G^(E) into G^ 2 \E) gives 



G^\x",x';E)=G ( - v \x",x';E) 
+ 



[G^ v \a 2 , a 2 ; E) - (c^(ai, a i; E) - ^ 

-G( v )(o 2 ,o 1 ;E)G( v )(o 1) o 2 ;E) 



x 



G^ v \ ai ,x';E) 



G^{x ,, ,a 2 ;E)G^ v \a 2 ,a 1 ;E) 



- ^ v \a 2 ,a 2 ;E) - ^^'(/.a^) 



+ G^ v \a 2 ,x';E) 



G iy) {x\a 1 ;E)G {y) {a u a 2 ;E) 



- (G<y\ ai ,a i; E) - ^jG( v \x",a 2 ;E) 



G^ v \x" , x'\ E) G<y\x", ai ;E) G< v > {x" , a 2 ; E) 

G( v \a u x';E) G^ v \ai, a±; E) — — G^ (a u a 2 ; E) 

7i 

G( v )(a 2 , z'; G^(a 2 , a i; £) G^(a 2 , a 2 ; £?) - — 

72 



1 



GW(oi,a i; £;)-l G^(a b a 2 ;E) 
7i 

G^)(a 2 ,a i; £) G^)(a 2 ,a 2 ;£)-- 

72 



(2.20) 



Bound states are determined by 



G^( ai , ai ;E n ) - — G^( ai ,a 2 ;E n ) 
7i 

G< v > (a 2 , ai; (a 2 , a 2 ; £ n ) - — 

72 



= 



The bound state wave-functions are given by 



K(x')*n(x") = lim 

Hi — >Hir, 



(E n -E)G^(x",x';E) 



(2.21) 



(2.22) 



Eq. (2.20) has an obvious generalization to the problem of N 5-functions. 
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We obtain 



G^ Sn \x",x';E) 

x(t")=x" 

dTe iET / h J Vx(t) 



1 

^ Jo 



exp 



x(t')=x' 



N 



m 



-x 



- V(x) + ^7fc<5(> - a k ) 



fe=i 





,x';E) GW(x", ai ;E) 


G( v \x",a N ;E) 

— ... G( v \a u a N ;E) 
7i 




,x';E) GW( ai , ai ;E)- 


G^(a N ,x';E) G( v \a N , ai ) 


... G^(a N ,a N ;E)- — 

In 




G( v \ ai , ai ;E)-- ... 

7i 


GW(a u a N ;E) 






GW(a N , ai ;E) 


GW(a N ,a N ;E)- — 
In 





dt 



(2.23) 



The quantization condition for the of bound states has the form 

G( v >(ai,ai;£ n )-— ... G^( ai ,a N ;E n ) 
7i 



G( v )(ojv,oi;f; n ) 
This can be written for short by 



G^ v \a N ,a N ;E n ) - 



In 



(2.24) 



det(ay) = , with the N x N matrix (ay) = (^G (1/) (a;, a,; £?) - y-J . 

' (2.25) 

The bound state wave-functions are given by 



K(x')*n(x") = lim 



(E n -E)G {SN \x",x';E) 



(2.26) 



Eq. (2.23) is proven in appendix A. Let us note that in the limit N — > oo Eq.(2.23) 
can be used to describe periodic 5-functions in solid state physics, see e.g. Goovaerts 
and Broeckx [29] for such a discussion, including the determination of energy bands. 
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2.3. 5 -Function Perturbation Along a Line 

The case of an arbitrary located two-dimensional 5-function perturbation on a line can 
now also be treated. We consider a 5-function perturbation oc — r yd(ax — by — c). This 
can be alternatively interpreted as two particles interacting pointwise with coupling 7. 
By a simple shift of variables and introducing center of mass and relative coordinates 
R and r, respectively 



r = ax — by — c , 

ab 2 x + a 2 (by + c) 



m 



a; 



R 



a 2 +b 2 



M = 771 



■+b 2 ' 

a 2 + b 2 
a 2 b 2 



(2.27) 
(2.28) 



we obtain by using the known exact propagator for the 5-function perturbation (e.g. 
[34] and references therein) the path integral identity 



x(t")=x" y(t")=y" 

Vx(t) J Vx(t)exp 

x(t')=x' V(t')=y' 
R(t")=R" 



= ab 



VR(t) exp 



R(t')=R' 
r(t )=r 

X 



iM 
~2hJt 



h J t i 



R 2 



m 



(x 2 + y 2 ) + r )5(ax — by — c) 



Vr(t) exp 



h L/ 



—f 2 + 7<5(r) 



dt 



r(t')=r' 



M 



2yri YiT 



exp 



M 



l\2 



x < exp 



2MT 

/J7 



(R" - R') 



^ 9 {r +r ) + -—kT 
h 2 h 2h 2 



+ 



2n 



dp exp 



dt 



p 2 hT 

i 

2/j 



x 



sinpr sinpr +cosjjr cos pr — 



e ip(r'+r") 

1 + i^ 

M7 



(2.29) 
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There seems no obvious way to incorporate more than one arbitrary line in IR 
equipped with 5-function perturbations. However, the two-dimensional result can 
be generalized to D dimensions. We introduce a 5-function perturbation on a hyper- 
plane a ■ x — 6 = (o,x£ M, D ). Introducing appropriate center of mass and relative 
coordinates [51, 52], we obtain similarly as in the two-dimensional case a (D — 1)- 
dimensional product of free particle kernels and one 5-function potential kernel in the 
relative coordinate. 
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2.4. 5 -Function Perturbations along Perpendicular Lines and Planes 

Let us consider first a two-dimensional system. Obviously we have the identity 



x(t")=x" 



y(t")=y' 



J Vx(t) J Vy(t)expl^J 



x(t')=x' 



y{t')=y' 



rn 



^(x 2 +y 2 )-V(x)-V(y) 



iVi 



N 2 



+ ^2 7i,fc*(E - a k) + ^2 ^^(y ~ hk ) 



fc=l 



k=l 



dt 



x(t")=x" 



i 



Vx(t) exp 



x(t')=x' 

y(t")=y" 



h .If 



m 



N i 



-X 



k=i 



dt 



x 



/ 



y(t')=y' 



Vy(t) exp { - j 



N 2 



k=l 



dt 



K^\x", x'; T) ■ (y", y':T) = K( N ^ (x", x', y", y'; T) . (2.30) 



The solution of each kernel in terms of its corresponding Green function is given by 
Eq. (2.23). Now call these Green functions G x Nl) (E) and G y N2) (E), respectively. The 
convolution theorem of the Fourier transformation now states that the Green function 
corresponding to the entire kernel G^ Nl ' N2 \x" , x\ y", y'; E) is given by the convolution 
of G x Nl) {E) and G y N2 \E), i.e. 



/oo 
dz G^ x Nl \x", x; z)G y N2 \y", y'\ E - z) . (2.31 
-oo 



) 



Obviously, we can generalize this result to the D-dimensional case. By separating the 
D path integrations we obtain a D-fold product of one-dimensional multiple 5-function 
potential kernels. The Green function of the entire kernel is given by a (D — l)-fold 
convolution of the corresponding one-dimensional multiple 5-function potential Green 
function, each of them given by Eq. (2.23). As complicated as this may look, for 
practical purposes only the case D = 3 is relevant, and what remains is a two- fold 
convolution. Because the number of 5-functions we take into account is not restricted, 
nor their variability in strength, we can take this model in the large Ni, N2 limit of 
perpendicular 5-function perturbations as a model for a rigid body lattice. The various 
entries in the determinants are simple, and all what one has to do, is the evaluation 
of the determinant and some numerical integration. Of course, the whole system can 
also be put into a (.D-dimensional) box (see next section) which only slightly increases 
the complexity of the system. A numerical investigation should yield a band structure 
of the energy levels. 
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2.5. Moving S -Function Perturbation 

Explicitly time-dependent problems in quantum mechanics turn out to be usually quite 
difficult and involved. In concrete physical situations one encounters mostly scattering 
problems and it is necessary to study the quantum mechanical properties of a system 
with a moving background, e.g. electrons scattered by moving ions. One useful model 
of such set-ups are 5-function perturbations, simulating point interactions (e.g. [49, 
64]). In this Subsection we study the quantum motion of a particle in the field of a 
moving 5-function perturbation [64]: 



V(x,t) = -~f5(x - vt) , 



(2.32) 



where v denotes the constant velocity of the 5-function. I present an alternative 
treatment as given by Duru [19] (compare also [36]) based directly on the summation 
of the corresponding perturbation expansion. The path integral now has the form 



K^(x",x';t",t') 

/,f/\ II 
X\t ) — X , 



rn 



x 2 + r )5{x — vt) 



dt 



x{t')=x' 



lim 



m 



n^oo V 27r i e?i 



N/2N-1 „oo f . N 

n / ^ ex p s 



3=1 



h 



rn 



— {xj - + e)8{xj - vtj) 



(xj = x(tj), tj = if + ej, with e = (t" - t')/N = T/N in the limit N 
perform the coordinate substitution 



y = x — vt , 



with 



x" - vt" 



y' = x' 



vt' . 



(2.33) 
oo). We first 

(2.34) 



This yields 
K^(x",x';t",t') 



/.lf\ ll 

y(* )=y 

I 

y(t')=y' 
m 

~ 2ihT 



exp 



i mv 2 T 

h 2 



m 



27rihT 



exp 



Vy(t) exp 



(x"-x') 2 



+ exp 



rn 



y 2 + mvy + r y8(y) 



2 

i mv 2 T 

h 2 



y f =x / — vt* 



K^(y\y'-t\t') 



y'=x' -vt' 
y" =x" — vt" 

(2.35) 
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where is given by 



K^{y'\y'-t\t') 

°° f\'y\ n r*" f t2 " 

=£(£) I *.■■•/, * 

x /!'<">(!,!, !/ ';i',i 1 )A-W(!/2.!/i;<2,ii)...A'<">(!,",!,„;i",<„) 



n=l 



n-1 



x K^^y'-t'M) • JJ ^(0,0;^,^-!) • tf<<V, 0; t", U 

J= 2 (2.36) 



Here K^ v \t",t') denotes the path integral 



y(t")=v" 



K^\y\y'-t\t')= J Vy(t)e W lj [^f + mvy^dt 



y{t')=v' 



m 



2n\hT 



exp 



(2.37) 



We find 



\y" — y'\ 1 — ~ — ^ . mv , „ 

W - l v /Z 2^E + 1 — (y" - y') 



(2.38) 



and we can apply the convolution theorem of the Fourier transformation yielding 



G^(y",y';E) = ^ 



exp 



\y'\ + \y 



+ i — (y" - y') 
n n 



7 m 



h\ 2 



(2.39) 



K(t",t') can be easily calculated by inverse Laplace transformation [21, Vol.1, p. 247] 



c- 1 



p -l/2( p l/2 + ^-l e -aVpj^ = e a/3+/3 2 t erfc ^_^ +/3 ^ ^ (240) 



yielding finally (re-inserting x' and x") 
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K {5 \x\x'-t\t') 



m 



2%ihT 



exp 



lm 
2KT 



(x" - x'f 



+ —2 exp < - 



2% 



h 



v(x — xt ) — v (x —vt)+m—T 



x exp 



^( K - Bn + K - ct 'l) + i^r 



h 



h 2h 



x erfc 



. ,x" - vt"\ + \x -vt'\ - -7T 
2 i hT V ' 'ft 



7777 

— o- exp 
ft 2 



lmv 



ft 



+ -L r^exp(-i^We-(") 



exp 



,1 11 ,,. , ... imv / „ , vT 

ip(\x" - vt"\ + \x' - vt'\) + — — \x" - x 

ft \ 2 



1 + i 



.phf 



7777 



We obtain one bound state 



(2.41) 



^ bound \x,t) = ^^exp 



777-7 , 1 TTIV , . 1 777 / 7 9 , 

& - ^1 + — £— - U*) + — I To + V )t 



h 



ft 



2% \h 2 



(2.42) 



which is the result of Ref. [64] , and the continuous states have the form 



^ cont -\x,t) = -Uexp 



2n 



ip 2 M 
2m 

exp 



x < e 



ipx 



2 1 



l mv i mv 

ip\x — vt\ H — (x — vt) H — — — — t 



1 + i 



. ph A 



7777 



> . (2.43) 



For 7j = the usual result is recovered. 
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3. Path Integrals for Boundary Problems 



3.1. Motion in Half-Spaces and Boxes 

In Eq. (2.6) we consider now the limit 7 — > —00 which has the effect that an 
impenetrable wall appears at x = a [13, 52]. We set lim 7 ^_ 00 G^(E) = G { ~ Wall \E) 
i.e. we obtain 



G^ Wall \x",x';E) 

=y °° dT ° iET/h 1 

= G^(x",x';E)- 



x(t')=x' 

G( v \x",a;E)G( v \a,x';E) 



in 



x 2 - V(x) 



dt 



G( v )(a,a;E) 
Bound states are determined by the equation 

G^ v \a,a;E n ) = . 

Again, the wave- functions are given by the residua of G(E) at E = E n , i.e. 



(3.1) 



(3.2) 



*„(x) = lim 



E — >E„ 



Eg — E 

G( v \a,a;E) 



1/2 



GW(x,a;E) 



(3.3) 



The Green functions G^ V \E) usually can be written as a product of two linearly 
independent solutions of the corresponding homogeneous Schrodinger equation, i.e. 
we can consider in general the Ansatz for G^ (E) 



G { y\x\x';E) = fEAsix^BEM , 



(3.4) 



where Je is a prefactor which depend only on E, and A^(x>) and Be{x<) are 
functions of the larger (smaller) of x',x", and depending, of course, on E as an 
parameter. We set the barrier at x = R and consider the motion for the (right-) 
half-space x > R. This yields 

G(*>«\x»,x';E) = fEA E { X> )B E { X< ) - (3.5) 



and the quantization condition has the form 



A En (R) = 
15 



(3.6) 



Similarly, if we put the barrier at x = a and consider the motion in the (left-) half-space 
x < a, we obtain 

&*«*VW;E) = f E A E (x > )B E (x < ) - f^ E (a)B E (x')B E (x") 

B E {a) 



and the quantization condition has the form 

BeM=0 . 



(3.8) 



Repeating the procedure for the double 5-function perturbation, we consider 
the limit lim^^-oo G {&2 \E) = G^ Box) {E) and obtain for the motion in the box 

a < x < b 



G^ Box \x",x'-E) 

= l d Te iET / h 
ft Jo 



/,/f\ it 
X\t ) —X 



V Bo xx{t) exp 



x(t')=x' 



h .If 



in 



x 2 - V(x) 



dt 



G< v >(x", 


x'-E) G (v \x",b;E) 


G( v \x" 


,a;E) 


G^ v \b,x';E) GW(b,b;E) 


G^ v \b, 


a-E) 


G {v \a,x'-E) G( v \a,b;E) 


G<y\a, 


a; E) 




G {v \b, b;E) G^ v \b, 


a-E) 








G^ v \a,b;E) G^ v \a, 


a- E) 







(3.9) 



The quantization condition for the (infinite number of) bound states has the form 

GW(b,b;E n ) G^(b,a;E n ) 
GW(a,b;E n ) G^(a,a;E n ) 

The bound state wave-functions are given by 

K(x')*n(x") = lim \(E n -E)G^ Box \x",x';E) 



(3.10) 



E — >E„ 



(3.11) 



Taking the same Ansatz (3.4) as before for the Green function G^ V '(E) we can simplify 
Eq. (3.9) into 

G<-VV; B ) = f E A E{l> )B E ( X<) + AE(a)B J E _ AE{V)BE(a) 
x I A E (x")B E (a) \B E (x')A E (b) - A E (x')B E (b)} 



+ B E (x")A E (b) A E (x')B E (a)-B E (x')A E (a) }, (3.12 
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and the quantization condition has the form 



A En (a)B En (b)-A En (b)B En (a) = 



For radial problems with a boundary at r = R we obtain similarly 



x(t )=x 



h 



dTe iET/H 



x(t')=x' 



T^Boxx{t) exp 



Y^G\ Box \r" ,r';E)S^{Q.")S^{Q.') 



h Jt> 



id 



x 2 - V(r) 



dt 



1=0 



where the radial Green function is given by 



G\ V) (R,R;E) 



The energy levels are determined by 



G^ V (R, R; E n ) = . 
For a radial ring with motion constraint by a < r < b we obtain 



G\ Rin9 \r",r; E) = 



G^\x",x'-E) G\ V \x",b;E) G[ V) (x",a;E) 
G[ V \b,x';E) G[ V \b,b;E) g[ V) (b, a; E) 
G[ V \a,x';E) g[ V) (a,b; E) G[ V \a,a;E) 



and the bound states energy levels are determined by 



G\ v \b,b;E) G{ V) (b,a;E) 
G[ V) (a,b;E) G[ V) (a,a;E) 



G\ V) (b,b;E n ) G[ V \b,a;E r , 



G { P{a,b-E n ) G\ v> (a,a;E n ) 
The bound state wave-functions can be obtained from 



= 



n n {x')^ ln {x") = lim 

Mi — >rj„ 



(E n -E)G[ Rin9 \x",x';E) 
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3.2. Boundaries Along Perpendicular Lines and Planes 

Because the Green functions for problems with Dirichlet boundary-conditions can be 
constructed from the corresponding Green functions for 5-function perturbations we 
can also discuss the case of multidimensional boxes. We consider Eq. (2.30). Since 
we are only interested in one domain which emerges form the limit 71/2, k ~^ ~ 00 ? it 
is sufficient to consider the case Ni = N 2 = 2. According to the convolution theorem 
this gives the Green function for a two-dimensional system in a box a\ < x < bi, 
a 2 <y <b 2 

- / dTe iET / h 
ft Jo 

/,//\ // / 4-" \ a 

x(t )—x y{t )=y 



X 



V B oxX{t) 



T>Boxy(t) exp 



x(t')=x' 



y(t')=y> 




in 



(x 2 +y 2 )-V(x)-V(y) 



dt 



/oo 
dzGi B °*\x\x'-z)G y B °*Xy\y'-E-z) 
-00 



(3.20) 



This results can also be generalized in an obvious way to higher dimensions with the 
restriction that no more perpendicular planes are allowed than twice the dimension 
of the space, i.e. ^(boundaries) < 2D. The entire Green function is then given by a 
(D-l)-fold convolution of Green functions (3.9). G^ V \E) itself can include 5-function 
perturbations, c.f. the remark at the end of the last section. 



3.3. Moving Boundaries 

Actually, the problem of a moving 5-function perturbation belongs to a larger class 
of explicitly time-dependent potentials. One possibility is to change a usual potential 
according to V(x) 1— > V(x/()/( 2 , where ( = ((t) = (at 2 + 2bt + c) 1 / 2 , and one derives 
the path integral identity [15, 37] 

x(t"y- 



Vx{t) exp 



m 



1 



■V 



x 



x(t')=x' 



= (C C) exp 



2* l 2 (t)' VCW 



2h [ x (" x c> 



dt 



^,v(^,^;r(0) , (3-21) 



with (' = ((f), (" = ((t"), etc., and we have set r{t" ) = f* dt/( 2 (t). Furthermore 
u' 2 = ac — b 2 and K^ 1 ,v denotes the path integral 

z(s")=z' 



K u , tV (z", z'; s") = J Vz{s) exp U J 



m . 2 m ,2 2 , r( \ 
~2 Z ~ ~2 U Z ~ W 



ds 



z(0)=z' 



(3.22) 
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Another class of time-dependent problems has a time-dependence according to V(x) 
V(x - f(t)). Here one gets [19] (q' = x' - /', f = f(t'), etc.) 



x(t")=x' 



i 



Vx(t) exp 



m 



-x' 



~ V(x - /(*)) 



dt 



x(t')=x' 



exp 




f"{x" - /") - f'{x' - /') + 



f 2 (t)dt 



K !v {q",q'-T) 



(3.23) 



with the path integral Kv V (T) given by 

T)n( t) PYn < 

h 



K fy (q",q';T) = J Vq(t)exp 

q(t')=q' 




in 



q 2 - V(q) - mf(t)q 



dt 



(3.24) 



Hence, in both cases the time- dependent problem can be rewritten in terms of a 
time-independent problem, however, additional terms appear in the Lagrangian. Con- 
sidering now time-dependent 5-function perturbations, the time dependence can be 
transformed away and one is left with the path integrals K u i j(r{t")) and Kj- S (T), 
respectively, with the 5-function perturbation located at the coordinate origin. Let us 
denote the corresponding Green functions by G^^siE) and Gv S (E). Considering the 
limiting process of making the strength of the 5-function perturbation infinite repul- 
sive, gives the corresponding Green functions G^i t waii(E) and G j Wall (E) according 
to Eq. (3.1), with the totally reflecting boundary located at the origin, which de- 
scribe still time-independent problems. The time- dependent problem can be obtained 
by applying Eqs. (3.21,3.23), provided we can explicitly evaluate K u i t wall( T (t")) and 
Kv WaU (T). It is obvious that this last step is a rather strong restriction which allows 
only for a limited number of explicitly soluble examples, say if we consider u' = and 
/ = 0, respectively. One then obtains [(MB) - moving boundary, D = 1] 



K( MB \x" ,x';t" ,t') 

= (c"cr 1/2 ex P 




TS~ I ' ' { ff \ 



(3.25) 



and similarly as for the case of Eq. (3.23). 



4. Examples 

In this section I am going to discuss several examples of the formalism presented in 
the previous section. I do not discuss examples with 5-function perturbations, this 
has been already done in Ref. [34]. 
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4.1. Boundary Problems in Half-Space 



4.1.1. Free Motion in Half-Space. As the first trivial example we consider for 
completeness the free motion in the half-space x > (IR + ). We have for the free 
motion in IR 



x(t")=x" 



I 



i m 



x(t')=x' 



'P.r(t) exp [ '-^j- j x z di 



m 



2n\hT 



exp 



m 



2h\T 



[x" -x'f 



(4.1) 



(4.2) 



Insertion of G^\E) into Eq. (3.1) gives the Green function of the motion in the half 
space IR+, which can be transformed into the time-dependent propagator yielding 



K^ + \x",x';T) 



m 



m , „ , 2 \ / m 



2-nihT 

which is the classical result, e.g. [25, 57]. 



exp '-2Tnr (x ~ x) r exp Ki^T 



(x"+x') 2 



(4.3) 



4.1.2. The Harmonic Oscillator. In order to apply our formalism we have to know the 
Green function of the harmonic oscillator (HO). It is not difficult to derive form the 
propagator 



x(t")=x" 



1 



Vx(t) exp 



x{t')=x' 



im 

2ft 



J (x 2 -oo 2 x 2 ^dt 



exp 



i muj 



(x' + x" ) coiuT - - 



27riftsinu;T 1 2% 
the corresponding Green function G^ HO \E), which has the form [4] 

G^ HO \x",x';E) = - 



sinu/T 



(4.4) 



m 



nTC'u V 2 hu 



x D i 



+— 



2mu 



2mu 



h ) 2 + r- \ \ h 

Here one can make use of the integral representation [30, p. 729], ai > a 2 



(4.5) 



JO 



coth 2i/ ^ exp 



x 



ai + a 2 



coshx 



l2i±{^/ai<x 2 s\rmx)dx 
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W^MM^M (4.6) 



and exploit some relations between the Whittaker functions W u ^(z) and M VjfJl (z), 
and the parabolic cylinder- functions D v (z), respectively. The Green functions of the 
harmonic oscillator in the half-spaces x > a and x < a are then given by Eq.(3.1) 
(compare also the case of the so-called double oscillator [4]), respectively Eqs. (3.5,3.7), 
and the bound state energy levels are determined by 



2 nxuj 

D_i + e il ( - v/— ^— « ) = for x < a , (4.7) 



2mu 



+ ( \/ —r~ a ) = for x > a . (4.8) 



4.1.3. Morse Potential and Liouville Quantum Mechanics. As the next example for 
motion in a half-space we consider the Morse potential (M) 

V ( M \x) = ^!(e 2 - -2a e*) (4.9) 

2777- 

(Vo > 0, a G IR constants). This potential has Nmo,x < olVq — -| bound states. The 
Green function can be explicitly calculated by the path integral formalism [8, 16, 31, 
55] and has the form 



x(t")=x" 



rn ±2 _Y^ 2x _ 2 
2 2m v ' 



dt 



h 

»(*')= 

m r(| + y/-2mE/%- aV ) 
~ V h 2 T(l + ^^E/Ti) e (-'+-")/2 ' ' avo.^f^ v ~ ' " " r "av ,^ 

The Green functions in the half-spaces x > a and x < a are then given by Eq. (3.1), 
respectively Eqs. (3.5,3.7), and the bound state energy levels are determined by 

M a v ,V=2^/H(2V e a ) = for x < a , (4.11) 
W olVo , V ^E: /h (2Voe a ) =0 for x > a . (4.12) 

Liouville (L) quantum mechanics can be studied by setting a = in V^ M \ We have 
[38] 

G^ L \x\x'-E) = ^I V ^ /h (V e*<)K v ^ E/h (V e*>) . (4.13) 

Again, the bound state energy-levels E n in the half-space x > a are determined by 

K^K: /h (Voe a )=0 . (4.14) 

Note that the function K- lv {z) for large values of v is an oscillating sin(z)-like function 
(e.g. [53], p.142). 
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4.1.4. The Wood Saxon Potential. We consider the Wood-Saxon potential. It is 
defined as [26] (6, R, V > constants) 



V(r) 



r > . 



1 +e (r-b)/R> 

The corresponding radial path integral formulation has the form (D = 3) 

r(t")=r" 



(4.15) 



Ki(r",r';T) 



— -2 4- ^0 

2 r l + e (r-6)/Ji 



r(t')=r' 



dt 



(4.16) 

This potential is used to describe the potential trough for the strong interaction near 
the nucleus for radial symmetric nuclei. There are only a finite number of bound states. 
Putting / = makes this path integral looking like the "smooth step" potential in IR 
with a barrier at r = 0. The path integral solution of the "smooth step" potential in 
the entire IR in turn is given by [33, 48] 



x(t")=x" 

dTe iET ' h J Vx(t) 

x(t')=x' 

r(mi)r(mi + l) 

h 2 r(mi + m 2 + l)r(mi -m 2 + 1) 



ft Jo 



exp 



m . 2 Vb 

2 X + 1 + e (x-b)/R 



dt 



2mR 



x 



x 



1 - tanh ?f=± 



1 — tanh 



m 1 -m 2 , m 1 +m 2 

1 + tanh -^-jjf 



2R 



1 + tanh 



2R 



m 1 +m,2 
2 



1 — tanh 

x 2-F1 mi, mi + l;mi - m 2 + 1; — 



1 + tanh 

x 2F1 mi, mi + l;mi + TO2 + 1; — 



(4.17) 



Here denote mi 2 = V2m R(V-E - V ± v r -E)/h. With a barrier at x = a, such 
that we consider motion in the half-space x > a, we obtain that the Green function 
of the Wood-Saxon potential is given by Eq. (3.1), and the bound state energy levels 
are determined by (with < \E n \ < Vq) 



1 — tanh a ~ b 

2 Fi(/? + iA,/? + iA+l;l + 2# 1=0. 



(4.18) 



Here denote (3 2 = -2mE n R 2 /h 2 , A 2 = 2m(E n + V )R 2 /h 2 . The claim of Duru [18] of 
having solved this problem cannot be seen as correct due to an improper handling of 
the boundary-conditions. 
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4.1.5. The Linear Potential. The problem of the linear potential (L) V(x) = kx in the 
half-space x > has a well-known result, i.e. the infinite number of bound state energy 
levels are determined by the zeros of the Airy- function [26]. The propagator of the 
linear potential in the entire IR is quite easy to obtain due to the fact that it belongs 
to the class of quadratic Lagrangians whose path integral solution are determined by 
the classical action. It is given by [25, 57] (x G IR) 



x(t")=x" 



I 



Vx(t) exp 



x(t')=x' 



m 



i /•*" /to .o , \ ' 
— / — x — kx \ dt 

hj t , \2 J 

i / to (x" — x') 2 



2n\TiT 



exp 



h V 2 



hT , . ... h T 
— (x' + x") - — — 
2 v ' 24m 



(4.19) 



The corresponding energy dependent Green function cannot be evaluated in an obvious 
direct way by means of a Fourier transformation from the propagator due to the 
rather nasty T 3 -dependence. However, as it is shown in Appendix B, a space-time 
transformation does the job. For the Green function we obtain 



G^\x\x>-E) = \^ 



n E 

X 'J 



nl/2 



X K 



1/3 



2mk 



x> - 



E 
~k 



3/2' 



h/3 



2mk 



x< - 



k 



(4.20) 



The Green functions in the half-spaces x > a and x < a are then given by Eq. (3.1), 
respectively Eqs. (3.5,3.7), and the bound state energy levels are determined by 
(n e INq, x > a) 



K 1/s 



2mk 



En 
~k 



3/2' 



Ai 



a 



E n \(2mk\ 1/3 ' 
k J U 2 J 



= 



(4.21) 



4.2. Boundary Problems in Boxes 



4.2.1. The Inhnite Potential Well. Again we start with a simple example, i.e. with 
the infinite potential well (IPW) . We need the Green function of the free particle, and 
insert it into Eq.(3.9). Performing the Fourier transformation with respect to E then 
gives the well-known result [57] (— b < x < b) 



TO 
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X 



\ ex p 



m 



2\KT 



(x" - x' + 4nbf 



1 

4b 



x" - x' 



4b 



— exp 
8m6 2 J 



m 



2\KT 



(x" + x' + 2(2n + l)bf 



-e. 



x" + x' 1 TtHT \ 



46 



+ 



2' 



8mb 2 



1 oo . 



n=l 



2 2 



7m 



cos -—a; cos 



8m6 2 J \2b 



2b 



-x 



(4.22) 



Here 6 3 denotes a Jacobi 6-function. 



4.2.2. The Harmonic Oscillator. This particular problem has attracted some attention 
by Aguilera-Navarro et al. [2], Consortini and Frieden [14], Marin and Cruz [54], and 
Vawter [62] . By means of the Green function for the harmonic oscillator we obtain for 
the Green function in the box a < x < b Eq. (3.15), and the bound state energy-levels 
are determined by 




(4.23) 



4.3. Boundary Problems in Radial Boxes 



4.3.1. The Radial Box. In order to treat the problem of the simple radial box we 
consider the radial free Green function which is given by 



r(i")=r" 

(r , r ,, )^^J°° dTe iET ^ n J Vr(t)exp 

r(t')=r' 



r(t")=r' 



3-D 1 



= (rV)^- / dTe iET / h 
* Jo 



if' /m , /(/ + P -2, 
h L V 2 2mr 2 



V l+ D~2{r 2 ]Vr(t)exp ( i^r J f r 2 dt 



r(t')=r' 



h 2 (r 



h 



(4.24) 



Therefore the Green function for the simple radial box < r < R is given by Eq. (3.15) 
and the energy levels are determined by 



(4.25) 
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i.e. by the zeros z\^ n of the modified Bessel function I\ (z). The corresponding problem 
with r > R describes the free motion in D dimensions with a impenetrable (hard) 
sphere at the origin. There are no bound states. 

This result allows us to describe the problem of the radial potential 

™ - (*■«) 

within, respectively outside, a radial box. The Green function of the problem without 
boundary is given by [it can be considered as an analytic continuation of Eq (4.24)]: 

G (^) (r // r / ;E) = |j KA ^ v ^^^ /A ^ v ^^^ ^ (427) 

and the Green function with boundary in the space < r < R is given by Eq. (3.15), 
and the bound state energy levels are determined by 

I x (y-2mE n j^ =0 . (4.28) 

This kind of problems appear e.g. by considering a sector of angel a which has a 
boundary at r = R, where A = Ire/a (/ G IN), i.e. we have motion inside a region 
which looks like a piece of a cake [11]. 

Vice versa, we obtain for the Green function for the potential V\ in the space 
r> R 



G ( r > R \r",r';E) = ^K x (^^ ^^x^/^mE 



R 



-^^ K ^_y^ r _ } . ( , 29) 

Hence, we describe the problem of the potential V\ outside a hard sphere (hard 
disc) of radius R. A model for such a Green function is the Green function for a 
Bohm-Aharonov solenoid outside a hard disc of radius R with flux $ inside. Then 
A = \v + q$/27xch\ {v E Z), with q the charge of the particle and c the velocity of 
light (e.g. [42, 43, 63] and references therein). 

4.3.2. The Radial Harmonic Oscillator. The model of the radial harmonic oscillator 
(RHO) can either be seen as a radial harmonic oscillator inside a radial box, or a radial 
harmonic oscillator with a hard sphere at the origin. We have for the path integral 
solution for the radial harmonic oscillator [17, 28, 56] 

r(t")=r" 

Vr(t) exp 



/ 



r(t')-- 



hL V 2 2mr 2 2 
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r{t")=r" 

~ I 

r(t')=r' 



ix l ! d-2 [r 2 ]Vr(t) exp 



2h 



<Adt 



exp 



^(r' 2 + r" 2 )cotcT 
i n 



, moor r 

J, , D-2 . 



(4.30) 



The corresponding Green function is given by [4, 5] (c.f. Eq. (4.6)) 
G\ RHO \r»y-E)- r[ ^ (/ + ^"^ )] 



^(rV") D / 2 r(/ + L>/2) 

/ mcu 2 \ / 

X yv e 1/,, D-aJ — rr- 7V I M B l f( , D-2J 
27iu7 ' 2 V ' 2 M fl ^ I 2Kuj ' 2 2 /\ 



mw 9 



(4.31) 



With a barrier at x = R we obtain the Green function of the radial harmonic oscillator 
for a radial box by Eq. (3.15), and the bound state energy levels are determined by 



Me_ 

2hu > 2 



for r < R , 



W E n 1 n | £>-2 
2Bl^ ' 2 



R A = for r > i? . 



(4.32) 
(4.33) 



4.3.3. The Coulomb Potential. We consider the Coulomb potential (C). The path 
integral reads (x G IR D ) 



x(t )=x 



K^ c \x'\x';T) 



Vx(t) exp 



— X 2 + — — 

2 \x\ 



dt 



(4.34) 



x(t')=x' 

Here only the Green function can be explicitly calculated. The radial Coulomb Green 
function has the form (see e.g. [4, 12, 35, 40, 41, 59] and references therein) 

Gf\r", r>- E) = (rV)-^^' + ^ + ^ 



h\ 2E (21 + D-2)! 



x W_ ip;Z+i¥ 2 



8m£ 



h 2 '>J M -^+^\]i % 

with p = (qiq2/h)^/m/2E. With a barrier at a; = R we obtain the Green function 
of the Coulomb potential for a radial box by Eq. (3.15), and the bound state energy 
levels are determined by 



SmE . 
— r K ) , (4.35) 



\d q 1 q 2 / m I \ P - 2 ( V ' — &TTlE n 



for r < R , 



(4.36) 



W^im > — —,,n-2 (J-8mE n ? J = for r > # . (4.37) 
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4.4. Boundary Problems in Radial Rings 

It is not difficult to generalize the preceeding results to the case of radial rings. The 
corresponding Green function in a radial ring a < r < R are given by Eq. (3.17) and 
the bound state energy levels are determined by: 

4.4.1. Simple Radial Ring. 




(4.38) 



4.4.2. Radial Harmonic Oscillator in a Ring. 

( 2\ W ( TD^ 

M E n 1 ,, , D-2 \ I — — a \W E n l( l4 _D-2-,\ 



W En i/7_i_D-2\( — r (X \ M E n 1 D _ 2 [-— R 2 ) . 4.39 



4.4.3. Coulomb Potential in a Ring. 




(4.40) 



4.5. Moving Boundary Problems 



4.5.1. The Infinite Potential Well. We consider the example of the infinite potential 
well with one boundary fixed at i = 0, and the other moving uniformly in time 
according to L(t) = L ((t) [15]. According to Eq. (3.21), the result then has for 
to' = the form 



K^ IPW \x",x'-t",t') = {CC P T — exp 



x 



e. 



2L 



im / „2(" , 2 C' 
— x x — 

2h I c" C 



2Lr 



2mLg 



e. 



2L r 



2mLl 



■ (4-41) 
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4.5.2. Moving Boundary in the Half-Space. We consider the half-space with the 
boundary moving uniformly in time according to Eq. (3.23) and we consider f(t) = vt. 
Hence / = 0, and for a free particle we obtain 



K^ MB \x" ,x';t" ,t') 



m 



2nihT exp I n 



im 



v(x" - vt") - v(x' - vt') + —T 



X 



exp 



im 

2nr' 



'-vT\ 2 )-e W (^\x" + x'-v(t" + t') 



11 , ./\|2 



(4.42) 



Generally, if the "mirror" principle can be applied, say for the harmonic oscillator, 
the result reads 



K^ MB \x",x';t",t')=e W {^ 



v{x" - vt") - v(x' - vt') + —T 



X 



K (x" -vt",x' -vt';T) - K (x" -vt",-(x' -vt');T)j . (4.43) 



5. Summary 

In this paper I have presented a comprehensive approach to the theory of boundary 
problems with Dirichlet boundary-conditions derived by an appropriate limiting 
process from 5-function perturbations in path integrals. 

First, I could derive from a path integral formulation and an exact summation 
of a perturbation expansion, closed formula? for the problem of multiple 5-function 
perturbations in one dimension. This result could be used to derive an integral 
representation for the problem of multiple perpendicular 5-function perturbations 
along lines, planes and hyperplanes in D dimensions, in the form of (D — l)-fold 
convolutions of the corresponding one-dimensional problem. 

I could also discuss the problems of a 5-function perturbation along a line, 
respectively its equivalent of two particles interacting via a 5-function, and a moving 
5-function perturbation. 

Second, I could consider boundary problems with Dirichlet boundary-conditions. 
They were derived form the ones with 5-function perturbations by letting the strength 
of the 5-function perturbations be infinite repulsive. This gave closed expressions 
for the one-dimensional motions in half-spaces, boxes, radial boxes, rings and mov- 
ing boundaries. The advanced problem of multidimensional perpendicular boundaries 
along perpendicular lines and planes was given similarly as for the 5-function perturba- 
tions, i.e. by convolutions of the corresponding Green functions of the one-dimensional 
case, however, only two boundaries are required for each dimension, thus simplifying 
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the corresponding one-dimensional Green functions, and the number of boundaries 
must not exceed (2 x dimension). 

These models of multiple 5-function perturbations and boundaries, respectively, 
along perpendicular lines and planes should serve as a numerically simple approach 
for the investigation of band structures. 

Several examples put our results in a better context and illustrated the formalism. 
Here some long unsolved problems in path integration could be treated successfully, for 
instance the Green function of the linear potential, and the Wood-Saxon potential in 
the half-space x > a. For the former, there are an infinite number of levels determined 
by the zeros of the Airy-function, whereas for the latter there are only a finite number 
of levels. For the linear potential it was required to calculate the Green function in 
the entire IR which was done in Appendix B. 

However, the question arises how to deal with other than Dirichlet boundary- 
conditions [9, 10, 13], which present only one possibility to define a self-adjoint 
extension of the (free) Hamiltonian with boundaries. It is known that an appropriate 
chosen 5-function perturbation can do the job; at least in case of the one-dimensional 
free particle in IR + the propagator with general boundary-conditions at x = can be 
explicitly computed [13]. It is derived from two perturbed constructively interfering 
kernels, i.e. a modified "mirror" principle is applied. As we know, this simple method 
does not work in the general case. All what we know in the general case is the 
corresponding Green function for the perturbed-, respectively boundary-problem. Of 
course, the same difficulties arise for the case of more than one boundary and in higher 
dimensions. These problems, however, will be discussed elsewhere. 
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Appendix 1. Green Function for Multiple 5-Function Perturbations 



In this appendix I want to show Eq. (2.23). The proof goes by induction. It is obvious 
that it is sufficient to show the validity of Eq. (2.23) for the denominator. Let us 
consider now the determinant (n G IN) det(a^) = 0, with the (n — 1) x (n — 1) matrix 



(a ij )=(G^(a i ,a f ,E)- 5 ^j . 



(Af.f) 



We know that Eq. (2.23) is valid for n = 1,2. We make the step form n — 1 to n. 
Doing this, every matrix entry is changed in the following way 



O'ij 1 0>ij — 



Therefore we obtain by expanding the determinant 
det (^|"~=i) 
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(A1.2) 
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The last terms with the power of the a nn factor less than n — 2 are all zero because 
two columns are always equal. The last non-vanishing term can be rearranged by a 
permutation where the first column is put into the last with increasing second index 
of the aij, etc. The emerging sum is now the definition of the expansion of a n x n 
matrix with the required feature, i.e. we have 

det (^IL"=i) = det ( a ^IL=i) ' ( AL4 ) 

which shows the induction. 



Appendix 2. Green Function for the Linear Potential 



We consider the path integral formulation for the linear potential 

x(t")=x" 

K{x\x'-T)= J Vx(t)exp 

x(t>)=x> 



h .If 



III 



— x 2 — kx)dt 
2 J 



(A2.1) 



We perform the combined space-time transformation (c.f. Refs. [39, 47, 60] and 
references therein) 



V=l*-f) 



3/2 



dt = -y~ 2 ^ds , 
9 



(A2.2) 



and introduce a new pseudo-time s" 

t" 

s " = lj t/ y 2/3 (t)dt, 



(A2.3) 



in this path integral, where the lattice implementation is given by Atj = ^y- 1/ ' 3 y J ^{ 3 
to guarantee a symmetric transformation with respect to initial and final coordinates. 
This gives the transformation formulae 



K(x",x';T) 



dE 
2tP\ 



iET/h n( n / 



G{x",x';E) 



G(x",x';E) = ^-(y'y'T 1/6 I ds" e~^" k ^ K(y\y'-s") , 



(A2.4) 



(A2.5) 



with the transformed path integral K(s") given by 



y(s")=y" 

K{y",y'-s")= j Vy(s)exp 

y(0)=y' 



h Jo 



m . 2 h 2 5 



8my 2 9 



ds 
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V(s")=y" 



i m 



Pi/3[V ]Vy(s)exp [ — / y ds 



y(0)=y' 



rrr, 771 



m 



2ifts 



(A2.6) 



Here use has been made of the radial path integral formulation (2.12) with the 
functional weight \i v together with the path integral solution (4.30) for inverse square 
radial potentials. Making use of the integral representation [30, p. 719] 



r e - a/x - bx Jv ^ dx = 2 j v J2a(Vb 2 + -b) 

Jo L 

we obtain for G(E) 



K„ 



2a(Vb 2 + c 2 + b) 



(A2.7) 



G(x",x';E) = ^ 
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E 
1 
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, (A2.8) 



(compare Steiner [60] for the case E = 0). The modified Bessel functions ii/ 3 and 
if 1/3 can be rewritten as Airy functions. One has to make use of the following relations 
[53, p.75] 



Ai(z) = 



Bi(z) 



1/2 



1/2 



1 (z 



7T V3 
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(A2.9) 
(A2.10) 



Then we obtain 



2\ 1/3 



h\hk 



Ai 



E\ (2mk 



x I Bi 



E\ (2mk 



x< - 



1/3 



v^Ai 



x< 



1/3- 



E 

T 



2mk\ 



1/3 



. (A2.ll) 



32 



References 

[1] S.Albeverio, F.Gesztesy, R. J.H0egh-Krohn and H.Holden, Solvable Models in Quantum Mechan- 
ics, Springer Verlag, Berlin, 1988 

[2] V.C.Aguilera-Navarro, H.Iwamoto, E.L.Koo and A.H.Zimerman, Quantum-Mechanical Solution 
for the Double Oscillator in a Box; Nuovo Cimento B 62 (1981) 91 

[3] J.-P.Antoine, F.Gesztesy and J.Shabani, Exactly Solvable Models of Sphere Interactions in 
Quantum Mechanics; J.Phys.A: Math. Gen. 20 (1987) 3687 

[4] V.L.Bakhrakh and S.I.Vetchinkin, Green's Functions of the Schrbdinger Equation for the Simplest 
Systems; Theor.Math.Phys. 6 (1971) 283 

[5] V.L.Bakhrakh, S.I.Vetchinkin and S.V.Khristenko, Green's Function of a Multidimensional 
Isotropic Harmonic Oscillator; Theor.Math.Phys. 12 (1972) 776 

[6] D.Bauch, The Path Integral for a Particle Moving in a 5-Function Potential; Nuovo Cimento B 
85 (1985) 118 

[7] M.Bordag, D.Henning and D.Robaschik, Quantum Field Theory With External Potentials 

Concentrated on Planes; J.Phys.A: Math.Gen. 25 (1992) 4483 
[8] P.Y.Cai, A.Inomata and R.Wilson, Path-Integral Treatment of the Morse Oscillator; Phys.Lett. 

A 96 (1983) 117 

[9] M.Carreau, The Functional Integral for a Free Particle on a Half-Plane; J.Math.Phys. 33 (1992) 
4139 

[10] M.Carreau, E.Farhi and S.Gutmann, Functional Integral for a Free Particle in a Box; Phys.Rev. 
D 42 (1990) 1194 

[11] L.Chetouani, A.Chouchaoui and T.F.Hammann, Path Integral Solution for a Particle Confined 

in a Region; J.Math.Phys. 31 (1990) 838 
[12] L.Chetouani and T.F.Hammann, Coulomb Green's Function, in an n-Dimensional Euclidean 

Space; J.Math.Phys. 27 (1986) 2944 
[13] T.E.Clark, R.Menikoff and D.H.Sharp, Quantum Mechanics on the Half-Line Using Path 

Integrals; Phys.Rev. D 22 (1980) 3012 
[14] A.Consortini and B.R.Frieden, Quantum-Mechanical Solution for the Simple Harmonic Oscillator 

in a Box; Nuovo Cimento B 35 (1976) 153 
[15] V.V.Dodonov, V.I.Man'ko and D.E.Nikonov, Exact Propagators for Time-Dependent Coulomb, 

Delta and Other Potentials; Phys.Lett. A 162 (1992) 359 
[16] I.H.Duru, Morse-Potential Green's Function With Path Integrals; Phys.Rev. D 28 (1983) 2689 
[17] I.H.Duru, On the Path Integral for the Potential V = ar~ 2 + br 2 ; Phys.Lett. A 112 (1985) 421 
[18] I.H.Duru, On the Path Integrations for the Wood-Saxon and Related Potentials; Phys.Lett. A 

119 (1986) 163 

[19] I.H.Duru, Quantum Treatment of a Class of Time-Dependent Potentials; J.Phys.A: Math.Gen. 
22 (1989) 4827 

[20] A.Erdelyi, W.Magnus, F.Oberhettinger and F.G.Tricomi (eds.), Higher Transcendental Functions 

Vol.11, McGraw Hill, New York, 1955 
[21] A.Erdelyi, W.Magnus, F.Oberhettinger and F.G.Tricomi (eds.), Tables of Integral Transforms, 

Vol.1, McGraw Hill, New York, 1954 
[22] L.A.Falkovsky, Density Attenuation of Surface Magnetic States; Sov.Phys.JETP 31 (1970) 981 
[23] L.A.Falkovsky and S.Klama, Broadening of Electron States in a Thin Film Within a Rough 

Surface in an External Longitudinal Magnetic Field; J.Phys.C: Solid State Phys. 20 (1987) 1751 



33 



[24] J.R.Farias, Simple Exact Solutions of the One-Dimensional Schrodinger Equation with Contin- 
uous Plus (5-Function Potentials of Arbitrary Position and Strength; Phys.Rev. A 22 (1980) 
765 

[25] R.P.Feynman and A.Hibbs, Quantum Mechanics and Path Integrals, McGraw Hill, New York, 
1965 

[26] S.Fliigge, Practical Quantum Mechanics, Vol.1; Die Grundlagen der mathematischen Wissen- 

schaften Vol.177 Springer- Verlag, Berlin, 1977 
[27] B.Gaveau and L.S.Schulman, Explicit Time-Dependent Schrodinger Propagators; J. Phys.A: 

Math.Gen. 19 (1986) 1833 
[28] M. J.Goovaerts, Path-Integral Evaluation of a Nonstationary Calogero Model; J.Math.Phys. 16 

(1975) 720 

[29] M. J.Goovaerts and F.Broeckx, Analytic Treatment of a Periodic 5-Function Potential in the Path 

Integral Formalism; SIAM J.Appl.Math. 45 (1985) 479 
[30] I.S.Gradshteyn and I.M.Ryzhik, Table of Integrals, Series, and Products, Academic Press, New 

York, 1980 

[31] C.Grosche, The Path Integral on the Poincare Upper Half-Plane With a Magnetic Field and for 

the Morse Potential; Ann.Phys.(N.Y.) 187 (1988) 110 
[32] C.Grosche, The Path Integral on the Poincare Disc, the Poincare Upper Half-Plane and on the 

Hyperbolic Strip; Fortschr.Phys. 38 (1990) 531 
[33] C.Grosche, Path Integral Solution of a Class of Potentials Related to the Poschl-Teller Potential; 

J. Phys.A: Math.Gen. 22 (1989) 5073 
[34] C.Grosche, Path Integrals for Potential Problems With 5-Function Perturbation; J. Phys.A: 

Math.Gen. 23 (1990) 5205 
[35] C.Grosche, Coulomb Potentials by Path-Integration; Fortschr.Phys. 40 (1992) 695 
[36] C.Grosche, Selberg Trace-Formulas in Mathematical Physics; Trieste preprint, SISSA/177/92/- 

FM, to appear in: Proceedings of the Workshop "From Classical to Quantum Chaos (1892- 

1992)", 21-24 July, 1992, Trieste, World Scientific, Singapore; eds.: G. Dell'Antonio, S. Fantoni 

and V. R. Manfredi 

[37] C.Grosche, Path Integral Solution of a Class of Explicitly Time-Dependent Potentials; Trieste 

preprint, SISSA/2/93/FM 
[38] C.Grosche and F.Steiner, The Path Integral on the Poincare Upper Half Plane and for Liouville 

Quantum Mechanics; Phys.Lett. A 123 (1987) 319 
[39] C.Grosche and F.Steiner, Path Integrals on Curved Manifolds; Zeitschr.Phys. C 36 (1987) 699 
[40] C.Grosche and F.Steiner, Classification of Solvable Feynman Path Integrals; DESY preprint 

DESY 92-189, to appear in the Proceedings of the "Fourth International Conference on Path 

Integrals from meV to MeV", May 1992, Tutzing, Germany, World Scientific, Singapore 
[41] C.Grosche and F.Steiner, Table of Feynman Path Integrals; to appear in: Springer Tracts in 

Modern Physics 

[42] A.Inomata, Remarks on the Experiment of Winding Number Dependence of the Aharonov-Bohm 

Effect; Phys.Lett. A 95 (1983) 176 
[43] A.Inomata and V.A.Singh, Path Integrals With a Periodic Constraint: Entangled Strings; 

J.Math.Phys. 19 (1978) 2318 
[44] A.Inomata and V.A.Singh, Path Integrals and Constraints: Particle in a Box; Phys.Lett. A 80 

(1980) 105 



34 



[45] W.Janke and H.Kleinert, Summing Paths for a Particle in a Box; Lett.Nuovo Cimento 25 (1979) 
297 

[46] S.Klama and U.Rossler, The Green's Function of Confined Electrons in an External Magnetic 
Field; Ann.Physik 1 (1992) 460 

[47] H.Kleinert, Path Integrals in Quantum Mechanics, Statistics and Polymer Physics, World 
Scientific, Singapore, 1990 

[48] H.Kleinert and I.Mustapic, Summing the Spectral Representations of Poschl- Teller and Rosen- 
Morse Fixed-Energy Amplitudes; J.Math.Phys. 33 (1992) 643 

[49] I.V.Komarov, Application of the Short Range Potential in the Calculations of the Ion-Ion 
Recombination; Sixth International Conference on the Physics of Electronics and Atomic 
Collisions, Massachusetts, 1969, abstract of papers, p. 1015 

[50] S.V.Lawande and K.V.Bhagwat, Feynman Propagator for the ^-Function Potential; Phys.Lett. 
A 131 (1988) 8 

[51] J.B.McGuire, Study of Exactly Soluble One-Dimensional iV-Body Problems; J.Math.Phys. 5 
(1964) 622 

[52] J.B.McGuire and C. A. Hurst, The Scattering of Three Impenetrable Particles in One Dimension; 

J.Math.Phys. 13 (1972) 1595 
[53] W.Magnus, F.Oberhettinger and R.P.Soni, Formulas and Theorems for the Special Functions of 

Mathematical Physics, Springer- Verlag, Berlin, 1966 
[54] J.L.Marin and S.A.Cruz, On the Harmonic Oscillator Inside an Infinite Potential Well; Amer.J. 

Phys. 56 (1988) 1134 

[55] K.Pak and I.Sokmen, A New Exact Path Integral Treatment of the Coulomb and the Morse 

Potential Problems; Phys.Lett. A 100 (1984) 327 
[56] D.Peak and A.Inomata, Summation Over Feynman Histories in Polar Coordinates; J.Math.Phys. 

10 (1969) 1422 

[57] L.S.Schulman, Techniques and Applications of Path Integration, John Wiley & Sons, New York, 
1981 

[58] M.Sieber, The Hyperbola Billiard: A Model for the Semiclassical Quantization of Chaotic 

Systems; DESY preprint, DESY 91-030 
[59] F.Steiner, Exact Path Integral Treatment of the Hydrogen Atom; Phys.Lett. A 106 (1984) 363 
[60] F.Steiner, Path Integrals in Polar Co-ordinates From eV to GeV; in "Bielefeld Encounters in 

Physics and Mathematics VII; Path Integrals From meV to MeV" , 1985, eds.: M. C. Gutzwiller 

et al. World Scientific, Singapore, 1986, p. 335 
[61] F.Steiner, On Selberg's Zeta Function for Compact Riemann Surfaces; Phys.Lett. 188B 

(1987), 447; Quantum Chaos and Geometry; in Recent Developments in Mathematical Physics, 

Conference Schladming 1987, eds.: H. Mitter, L. Pittner, Springer- Verlag, Berlin, 1987, p. 305 
[62] R.Vawter, Effects of Finite Boundaries on a One-Dimensional Harmonic Oscillator; Phys. Rev. 

174 (1968) 749 

[63] F.W.Wiegel, Path Integrals With Topological Constraints: Aharonov-Bohm Effect and Polymer 

Entanglements; Physica A 109 (1981) 609 
[64] S.K.Zhdanov and A.S.Chikhachev, Particle in a Field of Dispersing 5-Potentials; Sov.Phys.Dokl. 

19 (1975) 696 



35 



